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Abstract. We consider the method of infinite matrix inversion in the context of quantum 
state reconstruction. Using this method we give rigorous proofs for reconstruction formulas for 
. . . the Cahill-Glauber s-parametrized distributions and the rotated quadrature distributions. We 

also demonstrate how to construct the s-distributions from the quadrature data, 
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^ ■ 1. Introduction 

r-| '. Given a quantum system in a state p, that is, a positive operator of trace one acting on a 

Ph| Hilbert space Ti., one can perform measurements on the system to obtain probability measures 
^ • corresponding to various observables being measured. The inverse problem, namely, the recon- 
. struction of an unknown state from some set of measured probability distributions is one of the 
^ ! important problems in quantum theory, and consequently, it has been studied extensively. Here 
we are only interested in the case where the Hilbert space is infinite dimensional; this is typical 
in quantum optics (for an overview, see e.g. [42] or [24] and references therein). The "inversion" 
^ . character of the problem can be formulated in different mathematical forms. Probably the 
! most commonly used approach is to use some integral transform to convert the measured dis- 
^ I tributions into a desired quasiprobability distribution. The usual tomographic scheme uses the 
■ inverse Radon transform to reconstruct the Wigner function from the experimentally obtainable 
. rotated quadrature distributions [40l [H [38]. In the context of cavity QED and ion trapping, 
O I an alternative approach using an integral transform of the Rabi oscillations of a two-level atom 
coupled to the field has also been proposed [2H[26]. Even though the quasiprobability distribu- 
tions, such as the Wigner function, contain complete information about the quantum state of 
• r-j . the system, and can be used to calculate expectation values of observables, it is of interest to 
rN I also reconstruct directly the actual density operator of the system [TT], [25l [22], [TH [28l [Ml [32] . 
' In [25], an explicit reconstruction formula for the quadrature distributions was given. This was 
later generalized to cover the distributions of generic linearly transformed quadratures |13| . 
The recent progress in the field of quantum state reconstruction has been reviewed in |42) . 

The purpose of this paper is to illustrate the use of the method of infinite matrix inversion 
in reconstructing the density matrix of a quantum system from certain, measured or otherwise 
obtained, phase space distributions. This means that we are using a fixed computational basis of 
the (infinite dimensional) Hilbert space in question, and have to deal with various convergence 
issues related to the matrix inversion. These will then naturally lead to conditions on the state 
under which the formulas are valid. 

The structure of the paper is the following. In section[2]we fix the notations, describe the basic 
idea in the derivation of the reconstruction formulas, and discuss some general issues associated 
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with state reconstruction, infinite matrices, and approximating them by finite matrices. Section 
[3] is devoted to state reconstruction from the rotated quadrature distributions; we will give 
an alternative, mathematically rigorous derivation for the formula appearing in [25], and also 
demonstrate how only a finite number of quadratures need to be measured in the case where the 
density matrix is known to be finite. In section[4]we consider the Cahill-Glauber s-distributions. 
First we show how they can be constructed from quadrature data, and then we consider two 
methods for obtaining reconstruction formulas. Since the s-distributions are not in general 
positive, they are not strictly measurable quantities in the sense of the theory of measurement. 
However, they can indirectly be constructed from actual measurements. We return to this 
question in section 14.11 

2. Preliminaries on state reconstruction and infinite matrices 

Let H he a complex infinite-dimensional Hilbert space; we fix an orthonormal basis {\n) \ 
n e N} of TC. (Here we denote N := {0, 1, 2, . . .}.) This computational basis is identified with 
the photon number basis, or Fock basis, in the case where TC is associated with a single mode 
electromagnetic field. We will, without explicit indication, use the coordinate representation, 
in which H is represented as L^(M) via the unitary map H 3 \n) hn E L'^{^), where is 
the nth Hermite function. 

Let a and a* denote the raising and lowering operators associated with the above basis of 
and define the operators Q := -^{o* + a) and P := -^{a* — a) which, in the coordi- 
nate representation, are the usual multiplication and differentiation operators, respectively: 
{Qijj){x) = xip{x) and (P^)(x) = — i^(a;)Q In the case of the electromagnetic field, Q and P 
are called the quadrature amplitude operators of the field. The selfadjoint operator N := a*a is 
the (photon) number operator; it generates the phase shifting unitaries R{9) := e*^^, 9 G [0, 2n), 
and we can define the rotated quadrature operators Qe by 



We will also need the displacement operator D{z) = e^"* 2; G C, for which holds D{z)* = 
D{z)-^ = D{-z) and R{9)D{z)R{9)* = D{ze'^). The latter implies D{re'^) = R{9)D{r)R{9)\ 



which will be important. 

Let C{7{) be the set of bounded operators on 7Y, and TiTi) the set of trace class operators. 
The former space is equipped with the operator norm || ■ ||, and the latter with trace norm 
II • 111. Since Ti is associated with a quantum system, these sets have physical meaning: the 
states of the system are represented by positive operators p G T{l-L) with the unit trace. The 
state operator is uniquely determined by the numbers pnm '■= ('^IpI"^), constituting the density 
matrix in the number basis. The pure states correspond to projections onto the one-dimensional 
subspaces of TC, and are thus of the form \{p){ip\, where (p E TC is a unit vector. In particular, 
we will need the coherent states \z), defined for each 2; G C by 



To be precise, one has to take the operator closure when defining Q and P in terms of a and a*, but this is 
a well-known technical issue which is not important here. 
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Qe = R{9)QR{9y, 9e [0,27r). 
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When reconstructing an unknown state of the system directly from some measurement data, 
one of course requires that the data determines the state uniquely. The set of measured observ- 
ables is then said to be informationally complete [34l[6]. The observables in quantum mechanics 
are represented by normalized positive operator measures (POMs) E defined on a cr-algebra E 
of subsets of some outcome set The probability measure X ^-^ tr[pE(X)] corresponding 
to a state p and a POM E is then the one according to which the measurement outcomes are 
distributed. In this article, Q is always a subset of either R or = C, and S is the associated 
Borel cr-algebra B{Q). The most common observables are of the conventional von Neumann 
type, that is, normalized projection valued measures or, in the case when f2 = R, selfadjoint 
operators in 7i. In particular, the spectral measure of Qe is such an observable; we will denote 
it by Qe : -B(R) C{n). 

We consider the reconstruction of the state p by determining the elements of the density 
matrix (pnm)- Before proceeding to the derivation of the reconstruction formulas, we will make 
some general remarks on state reconstruction in infinite-dimensional spaces. 

A typical way of dealing with infinite matrices is to approximate them by finite matrices. In 
the case of a density operator p, a natural measure of approximation is the trace norm: for the 
projection Pp := ^^Ig P = 1) 2, 3, . . ., we can define 

for all states p G T(7i) such that tr(PppPp) = X]n=oPnn 0- Note that for any p G T(7i), 
there exists a smallest po ^ {1;2, . . .} such that ii{PppPp) ^ when p > Pq. The following 
well-known lemma implies that we can approximate a state p in the trace norm by a state pp, 
which has only finitely many non-zero matrix elements: 

Lemma 1. LetH be a separable Hilbert space, and {Pn)nm an increasing sequence of projections 
on H, converging strongly to the identity operator. Then for each trace class operator p on Ti, 
the sequences {PnP)ne'M, (p-Pn)neN, and (P„pP„)„gN converge to p in the trace norm. 

One can also use the trace of p to quantify how well a truncation pp approximates p: Sup- 
pose that we have determined the first p diagonal elements p„„, n = 0,1, ...,p — 1, and that 
Sn=o > 1 — e. Then, since tr p = 1 one must have p„„ < e for all n > p. The positivity 
of p implies |pmnP < PmmPnn for all n, m G N, so \pmn\ < e for all n > p or m > p. Obvi- 
ously, this has a simple physical meaning in the case of an electromagnetic field: the condition 
Yl^=pPnn < e just says that the probability of measuring a photon number larger than p in 
the state p is less than e. In many concrete applications, one can choose a finite p G N such 
that it is practically impossible (e ~ 0) to get a photon number larger than p; in any case, no 
measuring apparatus can detect arbitrary high energies. This then implies that pmn = for all 
n > p or m > p, suggesting that we may a priori assume that the state matrix is finite. 

However, typically the maximal detectable photon number is not actually known, or depends 
on the construction of the measuring apparatus; hence fixing it to a finite value represents an 
artificial truncation of a system which is genuinely infinite-dimensional. In view of the state 
reconstruction, this is particularly significant, since the state is not known beforehand. There 
is no way of knowing how large number p one must take so as to get Yl'^=p Pnn < e for a given 
e. 

^Normalized positive operator measure (POM) is a map E : S ^ ^0^) which is cr-additive in the weak 
operator topology, and has the property E(f7) — I (the identity operator). 
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Next we describe how to reduce the reconstruction problem into the infinite matrix inver- 
sion in the case of some phase space quasi-probabihty distributions. The distributions we are 
considering in this paper have densities of the form 



m,n=0 



where the fnm are some functions of the radial parameter r, and 6 is the angle coordinate. The 
general heuristic idea to obtain reconstruction formulas from such distributions is the following. 
Integrating against exponentials e^^^, for a fixed A; G N, with respect to 9 over [0, 27r), one is 
left with a single sum 

(1) 

W,,kir) := — / e'^'Wpir, e)de = V p^n / e^^'=+"-'"^Vn,n(r) dO = Vp„+fc,„/n,n+fc(r). 

2^ ^0 2^ n Jo n 

" m,n=0 " n=0 

We consider two methods of converting this to a infinite matrix equation: integrating ([T]) with 
respect to some functions gi^k of r to get 



(2) 



Pn+k,n / gi,k{f)fn,n+k{'r) dr, 

n=0 



or differentiating I times with respect to r at r = tq, which gives 
d^ °° rf' 

(3) -^i[9k{r)Wp^k{r)]\r=ro = Pn+fc,n^[^fc(r)/„,„+fc(r)]|r=ro, 

n=0 

where gk are some functions. With fixed k, both of these are now infinite matrix relations of 
the form 



oo 



(4) yi = y^^ainXn, 



n=0 



and the idea is to choose the functions gk,i or gk in such a way that this relation can be inverted 
as 



oo 



(5) = X] ^riiyi- 



=0 



Of course, one has to take care of the convergence of the series, and other technical details; this 
has to be done separately in each case. In particular, even if (ain) has a formal inverse (&nm), 
i.e. the relation 



n=0 n=0 
4 



holds, it is still not clear if (l5|) is true. In order to briefly demonstrate pathological situations 
that could in principle arise, we will give the following example at this point: Take 
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Clearly, these matrices are formal inverses of each other, and the relation (Jll) is well-deflned for 
any complex sequence The sequence {yi) is then given by yi = xi + xi+i. Consider the 

following two cases: 

(a) Let Xn = 1, n E N. Then ?/„ = 2, n G N, and the relation (JH]) does not make any 
sense, since the associated series does not converge for any n (its partial sums form the 
sequence 2, 0, 2, 0, 2, . . .). 

(b) Let Xn = (—1)"", n G N. Then ?/„ = 0, n G N, so the series in ([5]) is well-deflned, 
converging to for any n. But this is not equal to x„, so ([5]) does not hold. 

It is easy to see that in this simple example ^ holds exactly when lim„^oo a;„ = 0. For a 
detailed treatment of inflnite matrices, we refer the reader to [9]. 



3. Balanced homodyne detection and quadrature distributions 

We begin with a short review of the precise mathematical formulation of the balanced ho- 
modyne detection, which is a well-known scheme of measuring quadratures, and is used quite 
frequently in quantum optics (see e.g. [24]). 

The description of the homodyne measurement is not entirely straightforward: for any inter- 
val X C M, the probability tr[pQ0(X)] for the quadrature measurement can be obtained only as 
a limit of balanced homodyne detection measurements. The reader should consult [2lj for the 
basic description; we also wish to mention a mathematically rigorous proof given in |16| . 

Balanced homodyne detector consists of a beam splitter, with a pair of photon counters at 
the output ports. The signal light beam, and a coherent auxiliary beam are sent into the beam 
splitter, and the suitably scaled difference between the photon numbers at the output ports is 
recorded as an outcome of the measurement. The scale factor is the amplitude of the auxiliary 
beam. 

The idea in the mathematical description of [16] is the following. One considers a sequence 
of measurements, each with larger auxiliary held amplitude than the preceding one and such 
that the amplitudes grow without bound. From the data of each measurement, one calculates 
the probability that the outcome lies in X. The sequence of the probabilities thus generated will 
converge to tT[pQe{2)] in the limit where the amplitude tends to inflnity regardless of the state 
p. Hence, for any flxed state, and any flxed interval, the correct probability can be obtained 
with arbitrary high precision by using a sufflciently large amplitude. 

For each state p, and each flxed 9, the probability measure X i— > tr[pQ6i(X)] has a density 
function x i— > W^^{x, 6). It is easy to see that the function W^"^ : R x [0, 2ti) ^ C so deflned is 
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actually measurable, and Integra Now we can write 



(6) / g{x)W^^^{x)dx = ^pn+k,n / g{x)fn,n+k{x)dx, 

where 

:= ^ r ^''wf{x,e)de- 

^" Jo 

fnmijx) ■ hfii^X^hfYii^X^ , 

and (7 : M ^ C is any bounded measurable function. The technical details leading to (El) are 
straightforward (see e.g. [HI Lemma 5]). The relation (I6|) now corresponds to (JS]), once we 
have chosen suitable functions g = gi^k- 

By doing this, we will end up with the known reconstruction formula (see (fTTl) below) which, 
up to our knowledge, was first given by Leonhardt and D'Ariano [25], who wanted to replace 
the traditional and mathematically troublesome inverse Radon transform scheme with a direct 
reconstruction of the density matrix in terms of the quadrature data. The same formula was 
later obtained as a special case of more general group theoretical results by Cassinelli et al. [7] 

Instead of using their methods, we derive the formula by using matrix inverse relations. In 
our opinion, this method more explicitly illustrates the technique by which the matrix elements 
are picked out by the suitable averaging functions, namely the derivatives of the so called 
Dawson's integral. Most of the construction is already given in our recent paper [18], where 
we gave a direct proof for the fact that the set of quadrature observables is informationally 
complete. Using this method, we will also explicitly demonstrate the fact that if the state 
matrix is a priori assumed to contain only a finite number of nonzero elements, then it is only 
required to measure a finite number of quadratures. 

Now consider the function F : M — > M, defined simply hy Y [x) = 2-^ daw(x), where 

daw(x) = e"""'^ / e*^ dt 
Jo 

is the well-known Dawson's integral (see e.g. [I] pp. 298-299] or [39, Chapter 42]). We are 
interested in the derivatives Y^p\ p e N. Since they appear in many physical situations, such 
as in a series expansion of the so called Voigt function in spectroscopy (see e.g. jS^ p. 69-70]), 
the computational evaluation of y^''^ has been studied extensively (see e.g. [30 l [3l [29l [37| ). 

According to Lemma 1 of p8], each function Y^^^ is bounded. Hence we can use the functions 
gi,k = Since each W^j, is integrable, the integrals 

w^ir-= [ gi,ki^)w^iix) dx 

Jr 

are well-defined. Putting 

cin{k) := / gi,kix)hnix)hn+kix)dx= {n\Y^^+^^\Q)\n + k), 
Jr 

^That W^'^ may be chosen to be measurable can be seen by writing the density operator in the spectral 
representation and using the fact that 6 R{0) is weakly continuous. Integrability follows from Fubini's 
theorem. 
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we thus have the matrix relation 

I 

(7) W^l, = J2cin{k)Pn+k,n, 

n=0 

where the infinite sum is now reduced to a finite one because of the special properties of the 
Dawson's integral [HI Lemma 4]. The coefficients Q„(A;) can be evaluated analytically (see 
Appendix B for the calculation); the result is 



{n + k)\ \n J 



It is already clear that the elements Pn+k,n can be solved recursively from ([7]). In order to get 
an explicit formula, we notice that for any fixed A; G N, the matrix relation ([7]) now assumes 
the form 



(9) yi = 

where 



n=0 



n 



(10) y^ = (-iy;^777^, x^ = {-lt-^ 



2'(fc + /)! ' " ' ^ y (fc + n)! 

According to [36l P- 43], the relation ([9]) can be inverted to give 



Pn+k,n- 



which immediately yields an explicit reconstruction formula for the density matrix elements, in 

p,k,l ■ 



terms of the quantities W^'l , 



This is the same as formula (34) of Leonhardt et al. [25] . 

The finite case. The above reconstruction formula holds for all states p, so it is not necessary 
to assume it finite, in particular. However, when the state matrix is finite, the quadrature data 
required to determine it is naturally smaller; we close this section by considering the case of 
finite matrices. 

Let z G C and p be a positive integer. Since 



if follows that 



p, z = 1, 



t=o 
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when 2;^ = 1, z 7^ 1, that is, when z = e^'^'^ilP ^ g G Z \ pZ. For any t G Z, define 6'* := 2T[t/p. 
Now the points e^^p, t = 0,l,...,p— 1, divide the circle T into p parts and 

gGZ\pZ, 



Fix p and define for each A; G N and X G i3(R) an operator Vp(X) by 

^ i=0 

so that for a given state operator p, the experimentally obtainable density x 1— > W^^'^(x) of the 
measure X 1— > tr[pVp(X)] is of the form 



where W'^'^''^ is as before, so it involves only p different quadratures. Now 

1 

(nlV^mim) = (nlQmim)- Ve'("-'"+'')^^ X G 
and, hence, by defining s = n — m + k, 

(n|Q(X)|m), m = n + k mods, 



(n|V,^(X)|m) 



0, otherwise. 



Suppose now that pmn = when n, m > p. Then for any k G {0, ...,p — 1}, 

p— 1 p—l—k 

tr[pVj(X)] = J2 Pmn{n\\/liX)\m) = Pn+k,nHQiX)\n + k) . 



m,n=0 n=0 



Using the density VF^^'^, we get for any bounded measurable function g : M. C 

I g{x)hn{x)hn+k{x) dx. 

n=0 

This is in the same form as ([6]), except that the sum is already finite. Choosing again g = y('=+20^ 
/ G {0, . . . ,p — 1 — fc}, we thus get 

(12) iy;i^ = 5^Q„(fc)p„+,,„, 

n=0 

where 

Jr 

Comparing this with ([7j), we evidently get the reconstruction formula 

(13) Pn^k,n - (-1) y 1^ J ¥{kTlV.' 



where k E {0, . . . ,p — 1}, n G {0, 1, ... ,p — 1 — k}. Note that the quantities W^jfi only involve 
information from the quadratures P^^ corresponding to 6* G {6*^, . . . 0^'^}. 

4. Cahill-Glauber s-parametrized distributions 

The s-parametrized quasiprobability distributions were introduced in quantum optics as 
mathematical tools [HI [ID], but they have since become accessible also to direct measurements. A 
scheme for indirect determination of these distributions was suggested by Vogel and Risken [40] . 
and the pioneering experimental work was done by Smithey et al. [38]. The s-distributions also 
arise when considering realistic measurements, where the detectors are not assumed to be ideal 
[23l[12|. In that case, the measurement outcome statistics correspond to certain s-distributions, 
where the parameter s is related to the efficiency rj of the detectors, by s = 1 — 2/rj. For our 
purposes it is convenient to define a parameter A as A = This gives a bijective mapping 
on C \ {!}. We will use A as the parameter for these distributions throughout the paper, and 
for this reason we will call the s-parametrized quasiprobability distributions A-distributions. 

Let A G C, |A| < 1. Define a bounded operator i^'^ := (1 — A) XlfeLo norm 
= |1 — A|. If |A| < 1 then is a trace class operator with 

oo 

tr[K^] = (l-A)5^A'^ = l. 

fc=0 

For each A G C, |A| < 1, define a (weakly continuous bounded) function : [0, oo) x [0, 2n) — > 
C{H) by W^{r,6) = D{re^^)K^D{re^^)* , where D{re'^% is the displacement operator. For each 
state p define the A-distribution as the phase space distribution : [0, oo) x [0, 27r) C, 

oo 

W^{r,9) :=tr[pl^^(r,^)] = {I - \)^\\k\D{re'^y pD{re'^)\k). 

k=0 

Note that A = gives us the Q-function, and A = — 1 gives us the Wigner function (up to a 
constant scaling factor) of the state. 

If A > 0, then the operator density W'^ defines a covariant phase space observable 

B{C) 3 Q^\Z) :=- [ W\r,e)rdrd9 G C{n), 

TT Jz 

in which case the A-distribution can be measured via eight-port homodyne detection. The 
eight-port homodyne detector consists of two pairs of photon detectors, and the amplitude- 
scaled photon differences Di and D2 for each pair are recorded. See [2l| for the description 
of the setup measuring the Q-function, and note that in our case one has to use the state K''' 
instead of the vacuum input in one of the ports (we refer to [17] for a detailed description). For 
other values of A, the A-distribution is not a measurable quantity in the sense of the theory of 
measurement. However, we will later demonstrate that it can be constructed from the rotated 
quadrature distributions, which are obtainable via balanced homodyne detection. 

Let us consider now the operator valued function W'^. Using the fact that D{re^^) = 
R{6)D{r)R{6)* , a direct calculation gives us the matrix elements of W^{r,9) with respect 
to the number basis: 

{n\W\r,9)\m) = e'("-"^)^ir^„(r, ^), 

where 

K^Jr) := {n\D{r)K'D{rr\m). 
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By the formula of Cahill and Glauber [8j (for a detailed proof of the formula, see Appendix B), 
we get 



(14) K^.Jr) = W— (1 - ^)m-n+lg-(l-A)r2^m-n^n^rn-n(^^2 - A - X-y) 

V TDj' 

where 

is the associated Laguerre polynomial. Note that the function 

(15) a"l;^-" ((2 - A - X'^y) = a"l;^-" ( - (1 - A) VVa) = J2 ^ _ (1 - ^)'"^'" 

can also be defined at A = and the extension is smooth with respect to r and A. 

Since ||Vr^(r,^)|| < \K^\ = |1 - A| for all r e [0,oo), Q e [0,27r), and the mapping 6 
W'''{r,9) is weakly continuous, we can define the sesquilinear form 

(V^,^) ^ ^ j\'^\m\r,e)^) de, 

which is clearly bounded (with norm at most |1 — A|). Thus, for each A; G N, the operator 

W^{r) := — / e''''^W^{r,e)d9 
271" Jo 

is well defined as a weak integral. In addition, we have 

W^^,{r) := tr[pW,\r)] = ^ r e^'^'W^{r,9) d9, 

for each state p. Since 

oo oo 

W^{r,9)=tT[pW\r,9)]= J] pUn\W\r,9)\m) = ^ p^„e^("-")X'™(r), 



m,n=0 m,n=0 

we have 

'•27r 



(16) Wi^,{r) = ^Y1 P'-- / ^'^'^'""'^'KUr) d9 = J2 Pn+k,nKn+kir). 

-m,n=0 •'^ n=0 

This corresponds to equation jl]), and will be the starting point for the reconstruction of the 
state (see sections [43l and [44l below) . 

4.1. Constructing the A-distributions from quadrature data. In a recent paper [20] it 
was shown, that the Q-function of a state p can be constructed from the rotated quadrature 
distributions W^^ , by means of a generalized Markov kernel. That is, for each z = -^(g + zp), 
(g,p) G ]R^ there exists a function M^'^ : R x [0, 2ti) M, such that 

{z\p\z)= [ [ M^'P{x,9)Wf{x,9)dxd9, 



for all states p. In [33], the generalized Markov kernel for the A-distributions was constructed. 
We will briefly recall the results. 
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First, define Wp{r,9) =: Wp{z) =: Wp{q,p), where z = re*^ = -^{q + ip)- We are interested 
in finding a kernel Mf^, such that we can obtain the A-distributions from the quadrature 
distributions by integrating: 



/'27r n 

(17) W^iq,p)= / / Ml'^ix,e)Wf{x,9) 

Jo Jr 



dxdO 



where W^^ is the probabihty density related to the quadrature observable as before. It is 
sufficient to show the validity of the equation for coherent states p = \a){a\, a G C, in which 

case 



= {a\D{z)K'D{zr\a) = {a - z\K'\a - z) = {1 - >^'\{ 



a - z\-' 



fc=0 



(18) = (1 - A)e-I"-I^ V = (1 - A)e-"-"'-"^ 



^(l-A)|a-^|2 



k=0 



Now assume first that A G M, |A| < 1 and define the function M^'° : M x [0, 2-k) — > C by 



.0,0, 1 + /1 + A 



where Y{x) = 2£ daw(x) as before. Define the functions M^'^ : M x [0, 27r) C by 

Ml'P{x,e) := MI'°{x - qcosO - psme,e). 
The function M^'^ can be represented as a series of Hermite polynomials [33] 

(19) Mr (x, e) = {l-X)Y, S^(^^2.(x - a), 



k=0 



where a := q cos 6* + p sin = -^{q + ip)-, and the series converges absolutely. 

Since the absolute convergence of the series in equation (fT9l) does not depend on A being real, 



the function M^'^ can be defined for all A G C, |A| < 1. 
Defining u := -\/2 Re(ae~*^), we have 

/ W^'Ax,e)dx={a\Qe{X)\a) = {ae-''\Q{X)\ae~'') = ^ [ e-(^-")'rfx, 

for all X G i3(]R), which implies that 

" / Mr(x,^)<^„|(x,^)^ = ^ r / Mr(a:,^)e-(^-)^^. 
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Following the calculations in |20j, we get 



Jo Jr 27r 



2A; ! ^ V ^ / Jo 27r 

fc=0 ^ ^ «=0 ^ ^ 



(A-l)^|a-z|2'= 



k=0 

CO 

X)\a-z\^ 
^ 1 

k=0 



which shows that 



Thus, equation (fT7|l holds for all states 



4.2. Shifting the A-parameter. Suppose that for a given A G (—1, 1) one has obtained the 
A-distribution of some state p. One might be interested in finding the distributions for a 
different value, say A'. It turns out that in some cases this shifting of the A-parameter is needed 
for reconstructing the state. Indeed, the reconstruction formulas we derive in the next sections 
are only valid for certain values of A. Therefore, we present here the formulas for this shifting. 

First note that e L\R'^) for all states p and A G C, |A| < 1 (see e.g. [HI Lemma 3.1]). 
This implies that the Fourier transform of Wp, as well as its convolutions with other integrable 
functions, are well defined. Now let A G (—1,1). According to equation (fTSi) we have for a 
coherent state a G C 

W^^^^z) = {I - X)e-^'-'^\--r 
Putting a = -^{x + iy) and z = -^{q + ip), we get 

W^^.)H(g,l>) = (l-A)e"^[(^-)^-(-^)^l. 
Let A' G (-1, 1), A' > A, and define the function ^-a.a' : ^ C by 

9x,x'[q,P) 2n \'-\ 

12 



so that (?A,A' £ L^{M?) for all A, A'. Now a direct calculation shows that the distribution W^|a)(Q,| 
is the convolution of W^|a)(Q,| with the function Qx^y, that is 



iW\^){a\* 9x,x'){u,v) = / W^^^i^^^{u-q,v -p)gx,x'{q,p)dqdp 

1 (i_AO(i-A)^ r 

27r A' — A /Tn.2 



^1 _ y)e-^[("-")'+(^-")'l = W^'. 



|a>(a| 



for all (m, f) G M^. Since this is valid for all coherent states, we have = * gx,x', for all 
states p and A, A' G (—1, 1), A' > A. 

The A-distributions obtainable via measurements are those, for which A > 0, so it is of greater 
interest to find the inverse for the relation above. For this, we need to be able to invert the 
convolution transform, which of course puts some restrictions on the distributions in question. 



A'- 



Let A and A' be as before. Then, by the Fourier theory, we have W^' = * gx,x' = Stt W^-gx, 
The Fourier transform of gx,\' can easily be computed, and we get 

9x,x'{u,v) = —e 2(i-A){i-A)^ 

ZTT 

for all {u,v) G M^, showing that gx,x' is pointwise nonzero. Thus, we have = ^j-^- 



G L1(M2), we have 



1 [ ^ (u, v) 



tX' 
P 

(2vr)2 7iR2^ gx,x'{u,v) 
(20) = — I e*^«"+^"^e^(^^^^)(^^"'+"'^l^,^'(M,T;)dMdt; 

for almost all G M^. Note that the above integrability condition can, at least in principle, 
be tested directly on the measured distribution, and thus no a priori information on the state is 
needed. It is easy to see that there actually exist states for which the condition is satisfied. The 
significance of equation (l20ll is that if one measures W^' for some A' > 0, then the distributions 
for other values of A, may be calculated, provided that the integral exists. In particular, the 
Q-function and the distributions corresponding to the negative values of A can be obtained in 
this way. 

4.3. Reconstruction via integration. We will first derive a reconstruction formula which 
can be obtained by integrating the A-distribution with respect to suitable averaging functions. 
The formula can be seen as a direct consequence of the overlap relation for the A-distributions 
[8] (the overlap relation has also been treated in e.g. pp. 58-59]). A different form of the 
formula appears also in [28]. However, we will go through the mathematical details to ensure 
the validity of the formula. As it turns out, the formula only works for negative values of the 
parameter A. 

Let A G (—1,0), and for each n, m G N, let : [0, oo) — > C be the function defined by 

equation (fT4l) . Substituting A with A^^ in equation (fT4l) we see that also the functions 
are well defined for each n, m G N. We wish to integrate the function Wp f^{r)Kiij^f^{r) with 
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respect to r, and in order to do that we need to pay attention to some convergence issues. This 
is done in the following lemma (for the proof, see Appendix A). 

Lemma 2. Let p be a state and K^^_^i be the function defined in equation (fT4l) . Then 

POO r'^ 

/ K,kir)^titkir)r dr = Y, Pn+k,n / K;^.„+fc(r)ir,y;,(r)r dr 
Jo „=o 

for allkeN, \e (-1,0). 

Using the orthogonality relation 

" e-^x"L^(x)L^(a;) dx = ^[^ + ^ + 1] ^^^ 

nl 

for the associated Laguerre polynomials [ISl p. 844, 7.414(3)], we calculate 



I 



oo 



W^,,ir)K^,^,{r)rdr = ^ p.+... J , ^ ' , ^, (1 - A)'=+^(l - X^^^'X^X-^ 



X 



^2/Cg-(i-A)r2g-(i-A-i).2^fc((2 - A - A-i)r2)Lf ((2 - A - A-^)r2)r rfr 







2 g /^-^-V (n + Ji!(/ + ^)! '"" I -'^-^^nm^-) dx 



/ n\l\ ,„_, r[A; + n + l] ^ 

2 4^^"+''"V (« + ^)!(/ + A;)! n' 



n=0 



1 /! (/ + A;)! 

2^'+'=''(/ + A;)! /! 
1 

— ^Pl+k^l 

for each A;, / e N. This corresponds to equation ([2]), where the series on the right-hand side has 
been reduced to trivial. Hence we can immediately write the reconstruction formula for the 
matrix elements: 

/•oo 

(21) Pn+k,n = 2 iy;;,(r)ir^;+,(r)rrfr. 

Jo 

Note that this method works only for negative values of A, and therefore the state can not be 
directly reconstructed from measurement statistics. However, as we have shown in the previous 
sections, the required distributions can be constructed from measurement data, at least in some 
cases. 

To further illustrate the fact that the negativity of the A-parameter is in fact a necessary 
condition for the formula, we give the following counter example. Consider the vacuum state 
p = |0)(0|, for which we have 

W^|o)(o|(^^) = (l-A)e-(^-^)'■^ 
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If we attempt to reconstruct the only nonzero matrix element poo by directly using equation 
(l2T]l . we find that 



Poo = 2(1- A)(l- A-^) / e-(^-^)(i-^ ^'■ rdr 



oo 



which clearly diverges for A G (0, 1), since in that case (1 — A)(l — A~^) < 0. Thus, the formula 
works for all states only if A G (—1,0). In order to obtain a valid formula for other values of A, 
we need to use a different method. 

4.4. Reconstruction via difFerentiation. We now proceed to the reconstruction formula 
obtained via differentiation. We will start with a lemma, which is needed for inverting the 
infinite matrix identity (l26ll appearing in the derivation of the reconstruction formula. The 
proof can again be found in Appendix A. 

Lemma 3. Let {xn)nm be a sequence of complex numbers such that 



n=p 

(a) For each n eN, denote 

oo k 

" ^ ' ip - n)\in' - p)\ 

n'=k+l p=n ^ ^ ^' 



Then 



oo 



E 



(p — n)\ 

p=n ^ ' 

if and only if lim^^oo -Rfc = 0. 
(b) Suppose that lim„_,oo \xn\ = 0. Then limfc^oo -Rfc = for all n G N. 

With the previous notations we have 

oo 

(22) e(^-^)'-V;;,(r) = p.+.,„e(^-^>^ir„V/.(0, 

n=0 

where the series converges absolutely and uniformly on any finite interval [0,R). In order to 
be able to differentiate the series term by term, and evaluate the value at the origin, we must 
prove that the series for each derivative converges uniformly in an interval [0,e), where e > 0. 
This is done in the following lemma (for the proof, see Appendix A). 

Lemma 4. Let p be a state and A G C, |A| < 1. Then there exists an e > such that the series 



(23) E 



Pn+k,n' 



n=0 



converges uniformly on [0, e) for each I G N. 
The series in equation ( l22l ) can thus be differentiated termwise. We then get 



/ J Pn+k,n 
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dr^ 



5 

r=0 



which corresponds to equation (14.4p . Now 

Qlj,2u+k 



r=0 



SO that we get two (nonzero) cases: (i) when k is even, then / must be even, and (ii) when k is 
odd, then / must be odd. 

Let h, j eN and consider case (i) when k = 2h and 1 = 2]. Now 



Pn+2h,n\ 

n=max{0 J — /i} 



n\ (1 - A)2J+iA"--'+^(n + 2h)\{2])\ 
{n + 2h)\ {j - hy.{n - J + h)\{j + h)\ 



(24) = ( ,^M(l-A)2^+iA^-^- Pn+2MAV(^ + 2/^)'^! 



.J ^ - ux {n-j + h)\ 

n=max|0,_;— raj 

Consider then case (ii) and put k = 2h + 1 and I = 2j + 1. Then 

r9.^ WA _/^2j + lV ,x2j+23^/.~j V" Pn+2fe+l,nAV(n + 2/l+l)!n! 

(25) W^p,2.+i,2,+i-^._ J(l-A) A }^ (n-j + /^)! ' 

n=max{0,jr— h} 

Define, for all n, m G N, 

Then both equations (i24l) and f l25i ) reduce to the single one: 



rn+k,n Tjr\ / ^ 



where either k = 2h and I = 2j or k = 2h + 1 and / = 2j ' + 1. 

In order to invert equation (i26ll . we need to check that the necessary and sufficient condition 
of Lemma [3] is satisfied. It turns out that there are two different cases. If |A| < |, then the 
equation may be inverted regardless of the state in question. If A > |, then one may need some 
prior knowledge of the state, namely, that the condition of Lemma [S] is satisfied. 

Using the identity T.T=oi-'^f{l) = which holds for n > 1 [IF, p. 3, 0.151(4)], 

we get by straightforward calculation 

n'=k+l p=n ' ^ ^' 



in! — ri){ri! — k — \\\{k — n)\ 

n'=fc+l ^ ' ^ ' 



oo 



/F+TTfc)!(^7Tfc)! 



_^ (f + A; — n)(f — l)!(fc — n)! 
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Since \pv+i+k,v+k\ < 1 and n < k, we get the bound 



< 



lAI 



,^{v + i + k)\{v + ky. 



v=l 

k °° 



< 



{k-n)\ 
{k + l)\\\ 

(A; + /)!|A|'= 1 
{k-n)\ (1 - |A|)'+'=+i 

{k + l)\ 1 



v + l + k 

V 





A 




1 - 




A| 



where we have used the identity = Yl'^=o ("t")^"' ^ahd when |2|<1. If|A|<|, then 



|A| 
1-|A| 



(1-^) 

< 1, and thus Hmfc^oo -R^ = for all n G N. 



Now Lemma [3] may be used on equation (l26l ) to get 

f2p + k 



1^ = V 

-^n+k,n / J 



P 



(1 _ ^)-(2p+fe+l)^p_ 



Thus, we get the reconstruction formula 



Pn+k,' 



OO 



2p + k 



(27) 



n 



[n 



^)!^ W (2p + fc)!(l-A) 



(l-A)2p+'=+i(p-n)!"'''^'2p+'^ 



2p+fc+l 



Qy,2p+k 



r=0 



which holds for all n, G N, at least if |A| < |. A similar bound for the parameter has also 
appeared in previous works concerning state reconstruction from the A-distributions, see e.g. 
[231 [12]. 

In this case the reconstruction formula ( l27ll works also for the Q-function. In fact, defining 
Qp{z) := {z\p\z) = Wp{r,9), where z = re*^, we get 



(28 



Pn+k,'. 



{2n + k)\ ar2"+^' 



Jo 



r=0 



Remark 1. We close this section with a note on reconstructing the state from the Wigner 
function. In the above discussion we have not made any assumptions on the state of the 
system, and as a result, the values of the parameter A had to be restricted. In particular, 
neither of the methods gives us a reconstruction formula for the Wigner function, that is, for 
A = — 1. Since in both cases the restrictions are related to the convergence issues, it is clear 
that one can obtain a reconstruction formula for the Wigner function if one knows a priori 
that the matrix elements of the state fall off fast enough. For example, if one assumes that 

the series X^^o \Pn+k,n\ converges for all G N, the integration method gives a formula 
also for the Wigner function. Clearly the coherent state p = \z){z\ satisfies this condition. It is 
easy to check that the differentiation method also works for the Wigner function of a coherent 
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state. That is, the series (l23l) converges for A = — 1 and the condition of Lemma [3] is satisfied. 
Another assumption which allows the reconstruction of the state for all (nonzero) values of A, 
is the finiteness of the state matrix. In that case, the differentiation method works well. It is 
interesting to note that the method works regardless of the size of the state matrix, as long as 
it is finite. A similar observation can be made about the integration method in the case of the 
Wigner function. 



5. Conclusion 

We have considered the problem of state reconstruction in the cases where the measured dis- 
tributions are either rotated quadrature distributions or A-parametrized distributions. We have 
proved three reconstruction formulas using the method of infinite matrix inversion, while paying 
close attention to the mathematical details. In particular, we have analyzed the convergence 
issues which, in the case of the A-distributions, give restrictions concerning the values of A for 
which the reconstruction formulas are valid. For the quadrature distributions we have given an 
alternative proof for the known formula ( |TTi ). In the case of the A-distributions we derived two 
different formulas, namely equations f l2Ti ) and ( l27l ). obtained via two different methods. The 
integration method provides a formula, which works only for negative values of A, and therefore 
cannot be used directly on measurement statistics. However, as we have seen, there are ways 
to go around this problem. The A-distributions for the negative values of A can be constructed 
either from the rotated quadrature distributions by means of a generalized Markov kernel, or 
from another A-distribution, corresponding to a positive A, by inverting a certain convolution 
transform. The differentiation method has also its limitations with 1/2 as the upper bound for 
|A|. 
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Appendix A 

We have collected here the proofs of the three technical lemmas needed in section [H namely 
Lemmas [2l [3l and [H 

Proof of Lemma [21 Let k and define a; := (2 — A — A^^)r^. Define for each N eN the 
function . , : [0, oo) C by 



„k „—x T k ( ^\ T k 



Now the function x e ^L^{x)Lf{x) can be estimated [14], and we get 

|xV"L^(x)Lf(x)| < 2\n + 1) ■ ■ • (n + k)e-^\L^{x)\ = 2^ + ^^' g-f |Lf (x) 
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for all X G [0, oo). This implies 

\fkM\ 



2 V + 



N 



^ ^ Pn+k,n 1 



/! 



n=0 



2 V il + k)\ 



\M ' |Pn.+fc,n| 



(n + fc)! 



A'^a;'=e-^L^(x)Lf(x) 



|AnxV-L^(x)Lf(x)| 



/! 



n=0 



(ra + A;)! 



n+k,n I 



n! 



N 



< 2" 



k-l 



il + k)\ 



|A| '^|pn+Ml 



n=0 



> + A;)! 



e-^Lf{x) 



Since |p„+fc,n| < 1 for all n, /c G N, the series Yl'^=o \Pn+k,n\ IM'^'J ^"^f'" converges for all A; G N by 



the ratio test. That is, for each k eN, there exists a > such that Xl^Lo IPn+fc,™ 
Cfc, for all N eN. Thus we have 



(n+fc)! 



< 



l/i^,w(^)l<2 



fc-i 



A;! 



|A|-'C,e-t|Lf(x) 



(/ + A:)! 



for all G N. Since -C/f (x) is a polynomial of order /, the function e 2 |Lf (x)| is integrable over 
[0, 00), so the claim follows from the dominated convergence theorem. 

□ 



Proof of Lemma [H Fix n G N, and consider the partial sum 



{p — n)\ 



p=n 

Substitute the convergent series i/p = X]^=p m- )\ i'^to it: 



^ in — n'\\ ^ ( n' — 



k 00 



p=n n'=p ^ p=nn'=p \ / p=o n'=0 



n'=0 



,p=0 



n \ I p 
p J \n 



n J \p J n'\ 



where we have used the result Y^^=q{~^Y^^ (!^p^ (^) — ^n,n', which holds when k > n' [36l p. 4]. 
We have proved (a). To prove (b), we use a crude estimate 

00 k 

K\ < E E 



00 k 

< E E 



(p — n)!(n' — »)! (p — n)\(n' — k)\ 

n'=k+lp=n ^ ^ n'=k+lp=n \f J \ ' 



l-^n'+fc I 

n'=fc+l ^ ' n'=l 
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Since lim„^oo \xn\ = 0, the sequence (x„)„gN is bounded, i.e. there is a i^' > with < K 
for all n G N. Since ^^^^ ^ < cxd, we can use e.g. the dominated convergence theorem to get 



oo oo 



= y = y lim = lim y 

^ ^ k^oo n'\ k->oo ^ n'\ 

n'=l n'=l n'=l 

which proves that limfc^oo -Rfc = 0. The proof is complete. 

□ 



Proof of Lemma [H From eqs. ffT4l) and ffTSll we see that 



{n + A;)! 



u=0 



u\ \n — u 



(1-A) 



2u 



d r 



First, suppose that A 7^ 0. Since |A| < 1, we can choose positive numbers e and 5 such that 
(1 + < 1 and |(1 — A)^A~^|e^ < 5. Let Uq be the smallest natural number with 2uq >l — k. 
Then for m > mq we have 

Ql^2u+k 



{2u + k){2u + k-l)---{2u + k-l + l)r'^^+^~^ < (2m + A;)'e 



\l2u+k-l 



for all r G [0, e). Thus we have 



9r' 



< 



{n + ky. 



|l-A|'=+^|A|"(n + A;)! ^ 



1(1 - A)^A-^|"(2M + fc)' 
'u!(n — m)! (A; + m)! 



.2u+k-l 



< |l-A|'^+i|Ar(n + A;)!e^-'^ 



< ii-Al'^+ijAr^^^^le'^"' 



5" (2m + A;)' 
u\{n — u)\ {k + u)\ 

^rn\^{2u + ky 



u=Q 



u=Q 



uj {k + u)\ 



Now the function u ^— ^'^^^'^x is bounded by, say, M > 0, so 



<M|1-A|'=+^(1 + 5)"|A| 



{n + A;)! 



_k~l 



for all r G [0,e). Since (1 + 5)|A| < 1 the series ''"^f^' (1 + 5)"'|A|" converges by the ratio 

test. The claim follows from the fact that < 1 for all n, A; G N. 

Suppose now, that A = 0. Let hq be the smallest natural number with 2no > I — k. The 
series (l23l) reduces to 



1 



.U2n+k °° 



d r 



'n+k,n~ 



1 



V^!(n + A;)! dr^ ^^^m+M"^^,^^ ^ (2n + A; - / + 1)! 
which clearly converges on any interval [0, R), R > 0. The proof is complete. 



(2n + A;)! 



□ 
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Appendix B 



In this appendix, we prove two formulas which we used in the derivations of the reconstruction 
formulas. The first one gives the coefficients Q„(A;) needed in the quadrature reconstruction 
formula, and the second one is the formula of Cahill and Glauber, first proved in [8J. Since the 
original proof does not contain all the mathematical details, we reproduce it here somewhat 
more carefully. 

The function y : R — > M defined in section [3] is an analytic function, and the derivatives of 
Y have the series representations 



k=o ^ -r I- 



for all a; G M (see proof of Lemma 1 in the appendix of [H]). We use the following well-known 
result da p. 838, 7.375(1)]: for all m, n, / G N 



'■mnl 



^''^ m.n.t. _^ ^ _^ / ig even. 



Hrr,ix)Hn{x)Hi{x)e-''^dx = { {^,{m+n-l))\ {l{n+l-m))\ {l{l+m-n))l 

0, m + n + / is odd 



where the last case follows from the fact that HmHnHi is an odd function when m + n + / is 
odd. Note that Imni = (in the even case) also when I < \n — m\, since then either ^{n + l — m) 
or |(/ + m — n) is a negative integer. 

For any j G N the operator Y^^\Q) is bounded and we can calculate its matrix elements. 
Let ]? G N be fixed. For all m, n G N 



{m\Y'-^P\Q)\n) 
(m|F(2p+i)(Q)|^) 



- [{m+n)/2] / , Nfc/, , N, r 

2)P ^ {~lr{k + p)\Im,n,2k 



V7r2™+"m!n! ^ 2^(2/^)! 

/ o\j3 [(m+n-l)/2] / .Nfc/, , , iM r 

-(-2)P {-Ij^jk + p + 1)1 Irn,n,2k+1 

Vvr2™+"m!n! ^ 2*^(2fc + l)! 



where [c] means the integer part of c; e.g. [1.5] = 1, and [-0.5] = 0. Note that {m\Y^'^P\Q)\n) = 
when m + n is odd, and {m\Y^'^P^^\Q)\n) = when m + is even. As we will see, one gets 
even more zero elements. 

First, let m + n be even. Since {m\Y^^\Q)\n) = {n\Y^'^\Q)\m) we may assume that m < n 
and define u' = n — m > 0. Then n = m + u' and m + n = 2m + u' so that u' = 2u where 
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M G N. Hence, let m, u E N and calculate using [36l p. 8] 



m+u 



{m\Y^^P\Q)\m + 2u) = (-2)V^' (m + 2M)!^ ^ (fc + p)' 



{m — k + u)\ {k + -u)! (A; — u)\ 



-1)p+"2P(p + m)! 
-l)f+"+'"2P(p + M)!, 



,m -/)!(/ + 2m)!/! 



m! /m + M— p— 1 



(m + 2u)\ \ m 



ml I p — u 



[m + 2-u)! \ m 
= 0, if and only \i p>u and m > p — u. 

Secondly, let m + n be odd, m < n, and define u' = n — m > 0. Then n = m + u' and 
m + n = 2m + u' so that u' = 2u + 1 where u G N. Hence, let m, m G N and calculate 

(m|y(2p+i)(Q)|^ + 2M + 1) 

^ ^ ^ ^ ' f^^{m-k + uy.{k + u + l)\{k-uy. 

[ i) V^-l^ + ^« + iJ-Z^^^_^),(^^2M + l)!/! 



= i-ir--^r^'^\p + u + i)J, — , + - - ^ - 1^ 

= (_i)P+-+™+i2P+i/2(p + u + 1)! J ^^^^ ~ 

^ ^ ^ ]j {m + 2u + l)\ \ m J 

= 0, if and only if p > m and m > p — u. 
Now we know exactly the matrix elements of any Y^^\ 

We will now proceed to the proof of the formula of Cahill and Glauber, namely the relation 

^^r)(A*\m\ = p-(l-^)l^l'n - V\m~n+l-m-n^nrm-ni(. _ \\2| |2 



(29) {n\D{z)K''D{zy\m) = W— e"^^"^^!"! (1 - A)'"~"+^^'"-"A"L™-"(-(l - A)'|z|7A), 

V TDj. 

where n, m G N. For any A G C define 

:= (1-A)A^ = (l-A)^A"|n)( 

on the domain 



5^|Ar"Kn|y.)r<oo . 



n=0 
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It follows from the usual spectral theory that (K^)* = K^. In addition, 

(30) D{z)\k) e D{K^), keN,zeC. 

Indeed, we have 

\{n\D{z)\k)\' = ^e-^'^'\z\"^-'^\Lr\\z\')\', n>k, 
and \Ll-\\z\^)\' < (see e.g. [H p. 786, 22.14.13]), so that 

n=k ' n=k ''' ' n=0 ^ 

the last series converging by the ratio test. Because of (l30l) . the operator D{z)K^D{z)* is 
densely defined, with the domain containing all the number states. Hence, the left hand side 
of the relation fl29ll makes sense. 

By fl30ll . also the coherent states belong to the domain of D{z)K^D{z)* . Using fl30ll again, 
we can write 

oo 

{p\D{z)K^D{z)*\p) = {P\n){n\D{z)K^D{z)*\m){m\(3) 

m,n=0 

(31) = e-\^\' Yl ^y={n\D{z)K'Dizr\m). 

m,n=0 V 

To prove (l29ll . one first calculates 

oo 

{P\D{z)K'D{zr\P) = {P- z\K'\P -z) = (l-X)Y A"(/3 - z\n){n\P - z) 
= (1 A)c~^^'^~^''+'^~^''^ Y A" ~ f 



n=0 



\2n 



n=0 



(1 - A) exp[(A - 1)1^ - z\^] = exp[(A - + \z\^ - (3z - z^)] 

(1 - A) exp[-|/3p + (A - l)|^p] exp[/5((l - A)^ + A;^)] exp[(l - A)^;^] 

(1 - A) exp[-|/3p + (A - l)!^^] Y ^[(1 - A)^+ A^]'" exp[(l - \)z^ 



ml 

m=0 



(1 - A)e-l^l^+(^-i)l^l^ y ^[(1 - X)zr ( 1 + ^-^^exp \( iil^^M!^ 



(1 - A)z 



where one has to assume that z 7^ 0, A 7^ 0, and A 7^ 1. Now one can use the formula 
(1 + y)"^ exp[-xy] = ^^=0 ""(a^)?/", which holds for \y\ < I ^ p.l038, 8.975(2)], with 
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X = —iiLJ^ML and y = (izf)=- This then requires that |/?| < |(A ^ — l)z\. One obtains 

■((1-A)|^|)^\ 



l3\D(z)K^D(z)*\ 



m,n=0 

oo 



-A 



ni,n=0 v^-f^- Y ■ \ / 

Putting z = re*^, we get 



oo 

e-l/3P 



\^ e'("-m)(i_A)A/^e(^"^)'^'[(l-A)r]™-"A"L"*-" ( - 

±;^n \ V m! V 



m,n=0 

By comparing this with (ISTll . we get the result. 
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